Polariton excitation in epsilon-near-zero slabs: transient trapping of slow light 
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We numerically investigate the propagation of a spatially localized and quasi-monochromatic 
electromagnetic pulse through a slab with Lorentz dielectric response in the epsilon-near-zero regime, 
where the real part of the permittivity vanishes at the pulse carrier frequency. We show that the pulse 
is able to excite a set of virtual polariton modes supported by the slab, the excitation undergoing 
a generally slow damping due to absorption and radiation leakage. Our numerical and analytical 
approaches indicate that in its transient dynamics the electromagnetic field displays the very same 
enhancement of the field component perpendicular to the slab, as in the monochromatic regime. 
The transient trapping is inherently accompanied by a significantly reduced group velocity ensuing 
from the small dielectric permittivity, thus providing a novel platform for achieving control and 
manipulation of slow light. 



I. INTRODUCTION 

Physical mechanisms driving to slow and fast light 
have attracted considerable attention from the scien- 
tific community in the last decades [H-0|- The inher- 
ent interest in slow light comes from the long matter- 
radiation interaction time, which can lead to consider- 
able enhancement of all nonlinear processes that in turn 
may be exploited for active functionalities [3], e.g. all- 
optical switching and modulation Q. The nonlin- 
earity may also be enhanced by reducing the effective 
area in subwavelength silicon on insulator and plasmonic 
waveguides 0, H| , where tight confinement opens up pos- 
sibilities for miniaturized nonlinear applications [9|4ll|. 
Alternatively, extreme nonlinear dynamics [I2T fl3ll. en- 
hanced second and third harmonic generation |l4 , Il5j 
is predicted in epsilon-near-zero (ENZ) metamaterials, 
where the linear susceptibility is tailored in such a way 
that its modulus becomes comparable with the nonlinear 
counterpart. Boosting the nonlinearity of ENZ plasmonic 
channels can also lead to active control of tunneling [l6| , 
switching and bistable response [17] . ENZ metamateri- 
als have also been used for directive emission fl9l . 
cloaking [20j , energy squeezing in narrow channels |2l| 
and subwavelength imaging 

In all of the above mentioned mechanisms, the trade-off 
needed to achieve enhanced active functionalities is paid 
in terms of increased losses. As a result, the ENZ regime 
one usually invokes refers to the case where the real part 
of the susceptibility becomes very small, while its imag- 
inary part remains finite. Indeed, due to the stringent 
physical requirement of causality, Kramers-Kronig rela- 
tions impose that dispersion be inherently accompanied 
by loss and the dielectric susceptibility can not become 
rigorously null. The residual loss either limits or even 
prevents giant enhancement of coherent mechanisms, e.g. 



in second and third harmonic generation setups [Til . [l5| . 
Recently, in the context of surface plasmon polaritons, 
a method has been proposed to overcome the loss bar- 
rier for superlensing applications by loading the effect 
of loss into the time domain [24j . In our analytical cal- 
culations, we will use a similar approach to study the 
behavior of an electromagnetic pulse that scatters from 
a slab having a Lorentz dielectric response. Indeed, by 
considering non-monochromatic virtual modes with com- 
plex frequency [25[ , it is possible to drop off the effect of 
loss on the temporal dependence of the "mode" itself. In 
this complex frequency approach, it is possible to achieve 
the condition where the dielectric susceptibility exactly 
vanishes. Our formalism treats the dielectric polariza- 
tion of the medium as a generic Lorentz oscillator that, 
in the epsilon-equal-to-zero condition, encompasses lon- 
gitudinal collective oscillations of both electrons (volume 
plasmons) and ions (longitudinal phonons) that can not 
be excited by light [26| , Recently, the question whether 
or not volume plasmons can be excited or not by classical 
light has been revived [27Tj3l| . Some studies on Mie ex- 
tinction efficiencies reveal a maximum around the charac- 
teristic frequency where the dielectric susceptibility van- 
ishes, attributing the enhanced extinction to the excita- 
tion of volume plasmons [28], [3(| ■ Conversely, other sim- 
ilar studies identify the physical origin of the enhanced 
extinction in the excitation of leaky modes [27J, l29( . The 
latter interpretation is also supported by studies of the 
excitation of surface phonon polaritons in ENZ slabs [32h - 

In this manuscript we numerically investigate and ana- 
lytically interpret the scattering of a spatially and tempo- 
rally localized optical pulse from a dielectric slab in the 
ENZ regime. We used a finite difference time domain 
(FDTD) algorithm to solve the full vectorial Maxwell 
equations coupled to the Lorentz oscillator equation for 
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the dielectric polarization of the slab. We find that, if the 
carrier frequency of the optical pulse matches the ENZ 
condition, electromagnetic quasi-trapping occurs within 
the Lorentz slab since, after the pulse has passed through 
it, an elecromagnetic-polarization (polariton) oscillation 
persists and generally slowly damps out. We demon- 
strate that non-trivial ENZ features like the enhance- 
ment of the longitudinal electric field component are still 
observable in the time-domain. We also find that the 
above mentioned phenomenology is not observed for op- 
tical pulses with carrier frequencies far from the ENZ 
condition. Thus, in order to grasp the underpinning 
physical mechanisms responsible for transient trapping 
in the ENZ regime, we analytically investigate the scat- 
tering features of the Lorentz slab by studying the virtual 
leaky modes of the structure. We recognize that a set of 
polariton modes with reduced transverse group velocity 
(v g ~ c/100, where c is the speed of light in vacuum) 
is excited. Indeed, the plasma frequency plays the role 
of a cut-off frequency and polaritons in the ENZ regime 
are intrinsically characterized by a reduced group veloc- 
ity. Thus, we are able to interpret the transient trapping 
by means of the excitation of slow polariton modes that 
are damped off due to medium absorption and radiation 
leakage in the outer medium. 

The paper is organized as follows. In Sec. II we re- 
port the results of numerical finite difference time do- 
main (FDTD) simulations, comparing the distinct phe- 
nomenologies occuring in ENZ and standard dielectric 
regimes. In Sec. Ill we analytical investigate the virtual 
leaky modes of the structure, we address their properties 
and we discuss their role in the interpretation of numeri- 
cal results are developed in section III. In Sec. IV we draw 
our conclusions. 



II. FINITE DIFFERENCE TIME DOMAIN 
ANALYSIS OF THE TIME-DOMAIN ENZ 
REGIME 

A. Pulse scattering by a Lorentz slab 

Let us consider the scattering interaction sketched in 
Fig. 1(a), where an electromagnetic pulse is launched 
along the z-axis in vacuum and orthogonally impinges 
on the surface of a dielectric slab. The pulse is a Trans- 
verse Magnetic (TM) excitation, with electric E x (x, z,t), 
E z (x,z,t) and magnetic H y (x,z,t) field components. 
The initial profile profile of the transverse electric com- 
ponent is 



E x (x,z,t = 0) = E e 
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This field is spatially confined both along the x- and 
z- axis, w x and w z being its transverse and longitu- 
dinal widths, respectively, it is centered at the point 
(x, z) — (0, zq) and it is longitudinally modulated with 
period A = 2nc/co. Hereafter we will focus on very 





FIG. 1: (a) Interaction geometry of the pulse colliding onto 
the dielectric slab, (b) Real and imaginary part of the slab 
dielectric permittivity for the Lorentz parameters used in the 
FDTD analysis as a function of the wavelength. 



long pulses such that w z ^> c/{ui). Thus, the quasi- 
monochromatic condition Slo/uj <C 1 is satisfied, where 
Q is the pulse central frequency and Suj ~ c/w z is the 
spectral width. 

The dielectric slab has width L, it is centered at 
(x, z) — (0,0) and we assume that, in the presence of 
the external electric field E, the dynamics of its dielectric 
polarization is governed by the Lorentz oscillator model 
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where u e is the resonant angular frequency, 7 is the 
damping constant, f e is the oscillator strength and eo 
is the dielectric permittivity of vacuum. It is well 
known that Eq.([2]) leads, in the frequency domain, to 
the constitutive relation D w = eoe(w)E tlJ where f u = 
J^ 00 dte luJt f (t) is the Fourier transform of f(t), D — 
eoE + P is the displacement field vector and e = 1 + 
f e /(uj 2 — ijcj — uj 2 ) is the frequency dependent medium 
dielectric permittivity. The realistic model of Eq.([2]) is 
particularly accurate for describing the medium dielectric 
response to fields with frequencies close to the resonant 
frequency uj e (so that contributions due to other reso- 
nances can be neglected). Therefore, the Lorentz model 
is particularly suitable for our analysis since we are here 
concerned with quasi-monochromatic pulses whose car- 
rier frequency Co coincides with (or is close to) the fre- 
quency 
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(3) 

where the real part of the permittivity vanishes, i.e. 
Re [e(wo)] = 0, the so called epsilon-near-zero (ENZ) 
regime. 

Wc have performed the numerical analysis of the pulse- 
slab collision by means of a Finite Difference Time Do- 
main (FDTD) scheme where the polarization dynamics 
of Eq.([5]) are coupled to Maxwell equations for the TM 
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FIG. 2: Results of FDTD simulation pertaining the interaction of the pulse (with u> = Wo) with the slab. The fields E x , E z 
and the Fourier transform F[E X ] are captured at three different time steps corresponding to the three rows. In the first and 
the second time steps the pulse peak is just behind and beyond the slab, respectively, whereas in the third time step the pulse 
has completely left the slab. Note the large E z component within the slab (signature of the epsilon near zero regime) and the 
persisting and pulse-free polariton oscillation in the third time step. 



field. Specifically, in order to isolate the relevant phe- 
nomenology characterizing the ENZ regime we have ana- 
lyzed through FDTD simulations two different situations 
where the same dielectric slab is hit by two spatially equal 
pulses with different carrier frequencies Co: the first (pulse 
0) is such that Co — loq so that it is suitable to scan the 
slab behavior in its ENZ regime; the second (pulse 1) has 
Cu = uj\ for which Re[e(wi)] > so that it experiences 
standard dielectric behavior. 

In view of the generality and ubiquity of the Lorentz 
model of Eq.([2]), we have chosen for our numerical sim- 
ulations a realistic medium with Lorentz parameters 
oj e = 3.75 • 10 15 Hz, 7 = 1.50 • 10 12 sec" 1 and f e = 
2.25 • 10 30 sec -2 in order to deal with optical pulses in the 
visible spectrum, the resonant frequency corresponding 
to the wavelength A e = 2irc/uj e — 0.502 fim. For such 
parameters Eq.Q yields cjo = 4.03 • 10 15 Hz and we have 
set wi = 4.23 • 10 15 Hz (for which Re[e(o; 1 )] = 0.42), 
the two frequencies corresponding to the wavelengths 
Ao = 2ttc/loq = 0.466 yon and Ai = 2ttc/oji = 0.445 fim 
respectively. In Fig. 1(b) we plot the real and imaginary 
parts of the dielectric permittivity for the chosen Lorentz 
parameters as functions of the wavelength A = 2ttc/ui, 
indicating the carrier wavelengths Ao and Ai that char- 
acterize the two pulses. We have chosen a slab width 
L = 0.41 /jm to minimize the pulse propagation features 
and to effectively highlight the impact of the medium 
polarization on field dynamics. We set zq = —150 jim, 



so that the pulse peak reaches the slab 500 fs after it 
has been launched. Pulse widths are chosen so that 
w x = 1 fJ-m and w z — 60 /im: the former is compa- 
rable with the central wavelength of the pulse in order to 
provide the optical beam with a non negligible longitu- 
dinal field component E x (see below) whereas the latter 
corresponds to a temporal width T — w z /c — 200 fs 
and a spectral width Su> ~ 1/T = 5 • 10 12 Hz so that the 
pulses are in the quasi-monochromatic regime. 



B. Pulse scattering 

In Fig. 2 we report the main results of the FDTD sim- 
ulation dealing with the interaction of the pulse with 
carrier frequency Q = loq with the Lorentz slab repre- 
sented in the figure by the semitransparent rectangu- 
lar blocks. The first two columns of the figure con- 
tain the plots of E x {x, z,t), E z (x,z,t) as functions of 
{x, z) whereas the third contains the Fourier transform 
F[E x ](k x , z,t) — f_^dxe' lk * x E x (x,z,t) as a function of 
(k x ,z); each row of the figure corresponds to a selected 
simulation time step. At the first time step (first row of 
Fig. (2)), t — 371 fs, the pulse is fully interacting with 
the slab (the pulse peak being about to hit the slab at 
t = 500 fs) and the electromagnetic field is character- 
ized by standard reflection/transmission features; in par- 
ticular both the reflected and transmitted pulses have a 
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FIG. 3: Results of FDTD simulation pertaining the interaction of the pulse 1 (with ui = Wi) with the slab. The fields E x , E z 
and the Fourier transform F[E X ], for comparison purposes, are captured at same time steps considered in Fig. 2. Note that 
both E x and F[E X ] are bell-shaped, that magnitudes of E z within the slab and in vacuum are comparable and that no residual 
polariton oscillation has been produced by pulse passage. 



transverse bell-shaped spatial profile and accordingly the 
Fourier transform F[E X ] is peaked as well. Note however 
that, even in this early transient stage of the interac- 
tion, within the slab the longitudinal component E z is 
comparable with E x and much greater than its vacuum 
counterpart. Such an enhancement of the electric field 
component perpendicular to the slab is a feature typi- 
cally associated with monochromatic ENZ regime, aris- 
ing as a consequence of the continuity of the displacement 
field component perpendicular to the interface [3, EH ■ 
Therefore this is the first evidence that the ENZ regime 
can effectively be observed in thoroughly realistic Lorentz 
slabs by means of an equally realistic scattering interac- 
tion configuration. The second row of Fig. 2 corresponds 
to the time step t = 767 fs a time when the incoming 
pulse (if freely propagating) would have passed behind 
the slab (its temporal width being 200 fs). Note that 
the longitudinal component E z is even greater than the 
previous time step, testimony to the fact that the ENZ 
regime also occurs in the time-domain. The transverse 
component E x shows novel spatial features, even more 
evidently displayed by its Fourier transform F[E X ] which 
is no longer bell-shaped and characterized by a complex 
multi-structured profile. The third row of Fig. 2 considers 
a later time step t = 941 fs much longer than the time 
spent by the pulse to fully travel into the slab and leave it. 
At this time step, the longitudinal component E z is still 
very large within the slab and the transverse component 
E x displays novel and unexpected features: it is symmet- 



ric under the reflection z — > —z, it is not transversally 
bell-shaped and its Fourier transform F[E X ] displays two 
peaks at the sides of k x — 0. Such phenomenology can 
be interpreted only by assuming that the interaction of 
pulse with the slab is accompanied by the excitation 
of a polariton mode whose oscillation lasts a time much 
longer than the pulse-slab interaction time. 



C. Pulse f scattering 

In order to appreciate the novelty of the above dis- 
cussed time-domain ENZ phenomenology, we now dis- 
cuss the interaction of pulse 1 with the slab, its carrier 
frequency being associated to standard slab dielectric be- 
havior. In Fig. 3 we report the results of the FDTD sim- 
ulation relative to pulse 1 and, for comparison purposes, 
we have given Fig. 3 the same structure as Fig. 2 with the 
same fields at same time steps. Remarkably, both E x 
and F[E X ] are everywhere and always bell-shaped, while 
the magnitude of E z within the slab is comparable with 
its vacuum magnitude. At the last time step the slab 
hosts no residual polariton oscillation resulting from the 
pulse passage. This is precisely the standard expected 
phenomenology of the reflection and transmission of the 
pulse by a dielectric slab, and the comparison with the re- 
sults of Fig. 2 proves that the phenomenology it contains 
is a manifestation of the time-domain ENZ regime 
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FIG. 4: FDTD predictions about the fields E x and E z , for both pulse and pulse 1 as function of (z,t) at a fixed plane 
x — x p = 1.76 /i. Pulse triggers a novel mechanism of metastable light trapping since its passage produces a strong and 
damped polariton oscillation which is absent in the case of pulse 1. 



D. Transient trapping in the time-domain epsilon 
near zero regime 

In addition to the remarkable fact that the same fea- 
tures of the monochromatic ENZ regime characterize its 
time domain counterpart (e.g. the slab hosts a pro- 
nounced enhancement of the field E z ), the results dis- 
cussed in the previous sections also clearly reveal that 
the scattering situation leads to the unique excitation of 
a polariton mode. In order to show more explicitly such 
a phenomenology, in Fig. 4 we have plotted the fields E x 
and E z , for both pulse and pulse 1 as functions of (z, t) 
at a fixed plane x = x p = 1.76 [im. The evident fea- 
ture that emerges is that pulse (see subplots (a) and 
(b) of Fig. 4) produces a strong and damped electromag- 
netic self-oscillation persisting a time (about 5500 fs) 
much longer than the probing pulse duration (200 fs), 
self-oscillation which is conversely not produced by pulse 
1 (see subplots (c) and (d) of Fig. 4) , whose electromag- 
netic track fades within the slab just after it has left the 
medium (at about t = 1000 fs). We conclude that, in 
the ENZ regime, the pulse travelling through the slab 
triggers a novel mechanism of transient light trapping. 



III. THEORETICAL ANALYSIS OF 
TIME-DOMAIN ENZ REGIME 

A. Polariton virtual modes analysis 

From the above discussed phenomenology, it is evi- 
dent that a quasi-monochromatic pulse with a spectrum 
centered at the zero of the real part of the slab permit- 
tivity excites a polariton mode that lasts a time much 
longer than the pulse-slab interaction time. In order to 
rigorously prove this statement and gain deeper under- 
standing of the underpinning physical mechanisms that 



support the time-domain ENZ regime, in this section we 
analyze the exact quasi-steady modes (virtual modes) of 
the slab. In our analysis we fully take into account damp- 
ing processes, which include medium absorbtion and ra- 
diation leakage in vacuum, adopting the complex fre- 
quency approach (25|. We start our analysis from the 
curl Maxwell equations for TM fields 



dE z | 8E X 
dx 



dH„ 



dz 

dHy 

dz 
dx 



= -Mo- 
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dE* 
dt 
dE z 



dt 



dj\ 
dt 

m 

dt 



where the polarization P(x,z,t) — P x (x,z,t)e x + 
P z (x, z,t)e z satisfies Eq.([2]) within the slab (\z\ < L/2) 
and it vanishes outside the slab (|z| > L/2). We take the 
Ansatz Aj(x,z,t) = Re \aj(z)e % ( kxX ~^ lt >] for every field 
component (A = E, H, P, a = e, h,p and j — x, z), where 
k x is the (real) transverse wavevector and fi = to — iV is 
the complex angular frequency with T > so that only 
damping modes are considered. Owing to the mutual 
temporal evolution of the electromagnetic field (E, H) 
and of the polarization field P, the Ansatz effectively 
amounts to considering polariton virtual modes. The 
magnetic field can be expressed in terms of the electric 
field components h y = —(k x e z + id z e x ) / (fj^fi) so that 
Maxwell's equations reduce to 



de.. 
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e(n,z) 



dz 2 



Q2 

—o-e(Q, z) - k{ 



dz ' 



0. 



(5) 
(6) 



where e(Q,z) = e(fi)0(L/2 - \z\) + 9{\z\ - L/2) (9(z) 
being the Heaviside step function) is the z-dependent di- 
electric profile. It is worth stressing that the permittivity 
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e is evaluated at the complex frequency f2. The general 
solution of Eqs. (|5I6|) is explicitly given by 
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where K = y — k^, k ~ ^/^e(O) - fc2, C* is the arbi- 
trary mode amplitude, = ±1 is a parameter that dis- 
tinguishes the symmetry of the solutions and S = ±1 is 
another parameter selecting the sign of the exponentials 
in vacuum. By construction, the modal fields in Eqs.© 
already satisfy the continuity of the field component par- 
allel to the slab surface (e x ) at the interfaces x = ±L/2. 
The boundary conditions (BCs) for the continuity of the 
displacement field component (ee 2 )perpendicular to the 
interfaces x = ±L/2 yield the dispersion relation 



(Ke - Ek) e lkL = Q(Ke + Eh), 



(8) 



which provides the complex frequency f2 for every given 
slab thickness L and transverse wave vector k x . We have 
solved Eq.(|8|) numerically and obtained the allowed f2 
corresponding to different values of k x using the same 
slab thickness and Lorentz dispersive parameters of the 
slab considered in Section II. In Fig. 5 we plot the results 
for the case = — 1 in the complex plane f2 parametrized 
through the wavelength A and the damping constant T 
(i.e. f2 = 2-kc/X — iT), using circles and stars for the 
E = — 1 and E = 1 modes, respectively, and using the 
marker color to label the value of the corresponding k x . 
Note that the S = — 1 and S = 1 modes belongs to two 
different branches which are characterized by the fact 
that the S = — 1 modes have damping constant greater 
than the S = +1 modes. This property can be easily 
understood by considering the z-component of the non- 
oscillatory part of the Poynting vector S = E x H for 
z > L/2; 



(S z ) =S e - 2 [ EIm M-#)+ r ']I eo Re 



\C\ 2 . (9) 



For S = — 1, the energy outflows from the slab and the 
damping of the virtual mode is more rapid since it loses 
energy through both medium absorption and radiation 
leakage. Conversely, for S = 1 the electromagnetic en- 
ergy is dragged into the slab, thus partially compensating 
for the medium absorption and consequently decreasing 
the virtual mode damping time (note that there is also a 
point where r = on the S = 1 branch corresponding to 
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FIG. 5: Virtual modes of the slab considered in Sec. II for 
the symmetry O = —1 in the complex plane Q, = 2-kc/X — 
iT. Circles and stars label the S = — 1 and S = 1 modes 
whereas the marker color labels the corresponding k x value. 
The complex frequency Qp is such that t(Qp) = 0. The 
thin continuous line is the temporal spectrum of the incoming 
pulse reported in arbitrary units. 



the exact balance between medium absorption and radi- 
ation drag). It is remarkable that the S = — I and S = I 
branches intersect each other at point flp for k x — > 0. In 
this limit the dispersion relation of Eq.© (for = — I) 
reduces to 



(10) 



and is satisfied only if e(Slp) = 0. Starting from the 
Lorentz model, it is straightforward to prove that the per- 
mittivity vanishes at tip = y/e+^e^ 7 2 /4 — ^7/2 that, 
for the above used Lorentz dispersive parameters, yields 
A P = 27rc/Re(r2p) = 0.4667 [im and T P = -Im(fip) = 
8.3 ■ I0 -3 /s _1 , precisely matching the point of Fig. 5 
where the two branches E = ±1 intersect each other. In 
turn, the plasmonic flp at which the permittivity van- 
ishes plays a central role in the analysis of the virtual 
modes. For the complex frequencies f2 reported in Fig. 5, 
|e(f2)| < 0.06 and therefore, all the obtained modes with 
symmetry = — 1 imply the time-domain ENZ regime. 

In the same portion of the complex plane f2 we have 
numerically found no allowed modes for the symmetry 
= 1. This can be grasped by expanding both sides 
of Eq.© in Taylor series of e (since |e(fi)| -C 1); at the 
zcroth order we readily obtain e~ klL — —0 which is not 
consistent if = I (and which, on the other hand, yields 
k x = for = -1). 



B. Interpretation of FDTD results in terms of 
virtual polariton modes and slow-light regime 

Usually, within the standard real frequency approach, 
the solutions for the slab modes with E = — 1 are dis- 
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regarded since they are considered unphysical. Indeed, 
if \k x \ > oj/c, the solutions with H = +1 represent con- 
fined modes propagating along the ^-direction, while so- 
lutions with 3 = — 1 arc unbound modes that diverge 
at z — > ±oo. In addition, the introduction of the com- 
plex frequency introduces an inherent field singularity 
in the far past t — > — oo. Due to such intrinsic singu- 
larities, it is strictly impossible to rigorously excite a 
single virtual mode, its global existence on the whole 
space-time being unphysical. However, both singular- 
ities occur asymptotically and therefore virtual modes 
provide a very adequate description of the transient ENZ 
slab behavior occurring within a spatially bounded region 
and through a finite time lapse. In order to prove this 
statement and to basically provide a theoretical analyt- 
ical description of the transient light trapping discussed 
in Sec. 2D, in Fig. 5 we have superimposed the temporal 
spectrum profile of the incoming pulse of Eq.([TJ) (using 
the thin continuous line) on the complex-plane virtual 
modal structure. Note that, due to its wavelength band- 
width 5X = (27tc/ojq) 8uj — 5.8 • 10~ 4 fxm, the pulse spec- 
trum centered at Ao overlaps a limited portion of the 
considered complex frequency plane so that, specifically, 
the sole virtual modes with | A — Ao | < SX/2 are actu- 
ally excited by the considered pulse 0. From Fig. 5 it 
is evident that the excited virtual modes are character- 
ized by the transverse wavevector k x spanning the range 

1.1 [imT 1 < k x < 2.7 [imr 1 and that, due to the fi- 
nite bandwidth of the impinging pulse, the excited vir- 
tual modes with largest amplitude are those around the 
central transverse wavevector k x = 1.9 /^m, which cor- 
responds to A = 27rc/Re(r2) close to Ao = 0.466. This 
observation is in striking agreement with the results con- 
tained in Fig. 2, where one can see that the transverse 
Fourier transform of the field has, at the latest time step, 
two peaks centered at k x ~ 1.9 [im and k x ~ —1.9 \im 
whose width is of the order of 5k x ~ 0.8 fim. There- 
fore, when the pulse impinges onto the slab, it ex- 
cites precisely the virtual modes analytically predicted 
in Sec.IIIA, which are compatible with its spectral struc- 
ture. As a further validation of this statement, note that 
the virtual modes excited by the pulse ( | A — Ao | < 
SX/2) have a damping constant T spanning the range 

1.2 ■ lO^/s- 1 < T < 4.9 • lO^/s- 1 (see Fig.5), which 
corresponds to the extinction time r = 3/T spanning the 
range 609/s < r < 2604/s. Also this prediction based 
on the above virtual mode analysis is in striking agree- 
ment with the FDTD results since, by looking at panels 
(a) and (b) of Fig. 4, one can see that the electromag- 
netic excitation persists for a time of the order of 3000 fs, 
which is compatible with the maximum extinction time of 
the excited virtual modes. In addition, the spatial sym- 
metry of the = — 1 virtual polariton modes matches 
the numerical results displayed in Fig. 4: the transverse 
field component (e x ) is antisymmetric with respect to the 
z = axis, whereas the longitudinal field component (e z ) 
is symmetric. 

The final ingredient needed to thoroughly interpret the 



transient light trapping observed in FDTD simulations 
is related to the intrinsic slow-light nature of the phe- 
nomenon, which may be preliminarily grasped by con- 
sidering a bulk Lorentz medium, where transverse plane 
waves satisfy the dispersion relation k(uj) = (uj/c)^ e(w). 
Neglecting medium absorption (7 ~ 0), one finds that 
the phase velocity is v/ = oj/k and the group velocity 



dw/dfc is 



Vg(u) = C 



(ii) 



Thus, in the ENZ regime the phase velocity diverges 
Vf — > 00 whereas the group velocity tends to zero v g — > 0. 
Even though the subwavelength Lorentz slab used in our 
FDTD simulations is not a bulk medium and absorp- 
tion has not been neglected, the rough argument above 
still predicts the correct outcome. Indeed, by numerically 
solving the dispersion relation of Eq. (JSj) without neglect- 
ing losses one finds that, at the optical wavelength Ao, the 
transverse phase velocity of virtual polariton modes is su- 
perluminal Vf = u>/k x ~ 10c, while the transverse group 
velocity is extremely reduced v g — duj/dk x ~ c/100. 
For this reason, it is now clear how the virtual polari- 
ton modes, once excited, do not disperse quickly in the 
x-direction and remain quasi-trapped within the slab ow- 
ing to the tremendously reduced temporal dynamics. We 
conclude that the above described transient light trap- 
ping can be fully interpreted and physically understood 
by means of slow polariton modes supported by the slab. 



C. Volume plasmons 

Although the above discussed numerical and analyti- 
cal analysis of the transient light trapping characteriz- 
ing the time-domain ENZ regime is quite exhaustive, we 
now discuss its connection with the purely longitudinal 
modes, either volume plasmons (collective oscillations of 
electrons) or volume phonons (collective oscillations of 
ions), which the Lorentz medium can support. Hereafter 
we focus on volume plasmons, considering an unbounded 
bulk Lorentz medium where the TM electromagnetic and 
polarization dynamics are described by Eqs. (12141) . For the 
plane-wave Ansatz Aj(x,z,t) — Re [a J -e l ( fexa:+fe * 2: ~ ot )] , 
where A = E,H,P, a = e,h,p, j — x,z, k x ,k z are 
the (real) wavevector components and f2 is the generally 
complex frequency one gets 



[iQ^lhy — k x e z + k z e x , 
k z h y = Qe Q e(fl)e Xl 
-tte e(tt)e z , 



(12) 



k x hy 



where e(f2) is the dielectric permittivity with complex 
frequency Q. Volume plasmons are purely longitudinal 
electric oscillations owing to to the collective motion of 
electrons and are not accompanied by the generation of 
magnetic field, a feature that for plane waves amounts 
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to the collinearity of the wave vector k = k x e x + k z e z 
and the electric field E = k x ~E x + k z E z . Therefore, im- 
posing the condition k x E = 0, i.e. —k x e z + k z e x = 0, 
Eqs.O readily yield h y — and e(Q P ) = 0. Thus, vol- 
ume plasmons are inherently involved in the time-domain 
ENZ regime we are considering in this paper. However, 
it is worth noting that the virtual modes of the Lorentz 
slab are polaritons, entities fundamentally different from 
volume plasmons (or volume phonons). Indeed, a vol- 
ume plasmon is strictly characterized by the condition 
e(f2) = that implies the severe dispersion = fip and, 
in the presence of the slab boundaries at z — ±L/2, in- 
evitably leads to the inconsistency E z — > oo within the 
slab unless E z = in the outer medium. This is consis- 
tent with the well-known impossibility to excite volume 
plasmons by means of light. On the other hand, from 
Eq.([7]) one can see that the virtual polariton mode com- 
ponent E z neither vanishes outside the slab nor diverges 
within it. This is because for polaritons the dispersion 
relation of Eq.© is not as stiff as the volume plasmon 
dispersion CI = Qp and is satisfied also for ^= Qp. In 
addition the volume plasmon is a purely electric oscilla- 
tion with strictly null magnetic field, whereas the consid- 
ered virtual polariton modes are accompanied by a mag- 
netic field. In turn, even though TM polariton modes and 
volume plasmons occur in the same spectral region and 
are accidentally connected by the fact that in the limit 
L >> A the Lorentz slab is almost equivalent to a bulk 
medium, conceptually they are very distinct entities. In 
view of this, we conclude remarking that volume plas- 
mons can not be excited by classical light and that the 
absorption peak observed in experiments [HI, l3(ij is due 



to the excitation of virtual polariton modes, confirming 
the results given in Refs. [27l |29|. 



IV. CONCLUSIONS 

In conclusion we have investigated both numerically 
and analytically the properties of the time-domain ENZ 
regime. Specifically we have considered a dielectric slab 
whose polarization dynamics has been described through 
the realistic and ubiquitous Lorentz model and we have 
analyzed its interaction with quasi-monochromatic and 
spatially confined pulses with carrier frequencies close 
to the crossing point of the permittivity real part. The 
FDTD analysis has shown that the pulse is able to ex- 
cite a polarization-electromagnetic (polariton) oscillation 
which is damped and persists for a time generally longer 
than the effective time required by the pulse for passing 
through the slab. The underlying nature of this excita- 
tion has been elucidated through the analysis of the slab 
virtual modes that turn out to be located in a portion 
of the complex frequency plane close to the plasmonic 
frequency characterizing plasmon/phonons longitudinal 
volume excitations. Remarkably, due to this spectral 
property, both the group velocity and the transverse ve- 
locity (parallel to the slab) of each virtual mode turn out 
to be very small and therefore, the time-domain ENZ 
regime be naturally regarded as a novel platform for dis- 
cussing and investigating a plethora of slow-light phe- 
nomena. 
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